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We consider decay of the inflaton with a quartic potential coupled to other 

fields, including gravity, but restricted to spherical symmetry. We describe 
CN \ analytically an early, quasilinear regime, during which inflaton fluctuations 

and the metric functions are driven by nonlinear effects of the decay products. 

We present a detailed study of the leading nonlinear effects in this regime. 

Results of the quasilinear approximation, in its domain of applicability, are 
Tij- . found to be consistent with those of fully nonlinear lattice studies. We discuss 

^^ , how these results may be promoted to the full three dimensions. 
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I. INTRODUCTION 

Reheating is a crucial ingredient of inflationary cosmology [Q]. In most models of inflation, 
the quasi-exponential expansion leaves the universe very cold and empty, with energy stored 
in coherent oscillations of the inflaton field. Some kind of nonlinear process is required to 
convert that energy into energy of ordinary particles. 

Products of the inflaton decay (which will be denoted as x but in general will include 
the short-wave fluctuations of the inflaton field itself) can cause metric perturbations in a 
variety of ways. Two possible processes are represented by diagrams of Fig. |l]. In the first 
process, x fluctuations scatter off the inflaton's homogeneous component (which is still large 
at that time), thereby producing inflaton fluctuations, which in turn cause perturbations of 
the metric. In the second process, x fluctuations drive the metric perturbations directly. 

We want to stress that the processes shown in the diagrams are quite generic and exist 
regardless of the particular mechanism of reheating. They exist even if the inflaton only 
decays perturbatively, although the scattering process of Fig. la is absent from the early- 
days perturbative treatment 0, in which the decay is represented by an effective friction 
term in the equation of motion for the inflaton field 0. Fig. la corresponds to a driving force 
in the equation of motion for that may allow fluctuations of to reach values Sip ~ 0. 
The question remains of course (which we will address in this paper) what portion of this 
large fluctuation resides in super-Hubble modes, in particular those relevant to structure 
formation. 

Recent progress in the theory of reheating came after recognition that in many infla- 
tionary models the early stages of reheating can be studied rather completely. These are 
models where reheating starts by a parametric resonance 0,^. Parametrically amplified 
fluctuations are almost classical, and nonlinear effects can be explored numerically in the 
classical approximation [Q. It was found that in many cases the resonant stage is followed 
by a stage (called semiclassical thermalization), during which nonlinear interactions of fluc- 
tuations become increasingly important, and power spectra of the fields eventually become 
smooth 1^-0]. Importance of the scattering process of Fig. la was recognized early on in 
lattice studies in rigid spacetimes ^, and since then the estimate 50 ~ (0.1-1)0 has been 
confirmed for a number of models |17|-[T0| . 



Meanwhile, applications of the linear theory of metric perturbations to a special type 
of models have shown that in principle also super-Hubble modes, which necessarily acti- 
vate gravity, can undergo growth at reheating [pHl^. It is therefore became important to 
combine the two lines of research and see if a similar super-Hubble growth can occur due 
to nonlinear effects. The first attempts is that direction were made by studies in restricted 
geometries: planar geometry in Ref . [|T^ and spherical symmetry in Ref. |T^ . Both studies 
have indeed found such a growth [|l^,0 . 



The main purpose of the present paper is to supply an analytical counterpart to the 
numerical method of Ref. [|1^]. Although our approximate method, which we call quasilinear 
theory, cannot describe the fully nonlinear stage of the evolution, it allows us to describe 
the relevant nonlinear effects at earlier times, when the nonlinearity is still weak. This 
method retains linear terms in fluctuations of and the metric and quadratic terms in x- 
It generalizes the method developed in Refs. 0,|l3] for rigid spacetimes to the case when 
metric perturbations are present. In particular, it will allow us to take into account the 



processes shown in Fig. 1, which are not included in the standard hnear theory of metric 
perturbations. (In this paper, we consider the case when the expectation value of x is 
initially zero, so these nonlinear effects are the leading effects of x-) 

The system of spherically-symmetric equations is written down in Sect. |I|. In Sect. 
pni , we construct the quasilinear theory, restricting our attention to a model with a quartic 
infiaton potential. Within the quasilinear theory, we identify two potential sources of metric 
perturbations at super-Hubble scales. The first is the "erosion" of the infiaton by x fluctua- 
tions, as shown in Fig. la; the second, is the metric perturbations caused by the anisotropic 
stress, a purely relativistic effect that can be viewed as a part of Fig. lb. In the subsequent 
sections we analyze these effects in turn. We find that unless the decay products themselves 
have infrared-singular spectra the "erosion" is not infrared singular and thus cannot be im- 
portant on scales relevant to structure formation. On the other hand, the anisotropic stress 
gives rise to long-range (secular) terms in the gravitational potentials and is important at 
large scales, in a calculation restricted to spherical symmetry. In the concluding Sect. [Vllij, 



we discuss the overall physical picture that emerges when we try to promote our findings to 
the full three dimensions. 



II. GRAVITATING SCALAR FIELDS IN SPHERICAL SYMMETRY 

Early on in reheating, the particle density is small, and a weakly interacting particle 
species has a very large mean-free path. For such a species, a spherically-symmetric field 
will be relatively stable with respect to fragmentation into higher partial waves. This means 
that the restriction to spherical symmetry, when it is desired or unavoidable, can be imposed 
directly on the fields (and the metric), rather than on compound quantities, such as the stress 
tensor. In this sense, our calculations will be quite different from those with spherically- 
symmetric fluids: in the latter case particles move in all directions; it is only the large-scale 
motion of the fluid that is radial. In our case, because of the large mean-free paths, the fluid 
approximation is not readily available. 

We consider a system of scalar fields minimally coupled to gravity and described by the 
Lagrangian 



S = d'^Xy/^ 



^-lT.9,udy.d''<P,-V 



(1) 



where i labels different fields, k^ = 8ttG, and V = V{(f)i) is the potential which can depend 
on all of them. The most general spherically symmetric metric is 

ds"^ = {-N^ + ilj^(3^)dt^ + 2^^f3dtdr + ip^dr^ + a^r^dQ^ (2) 

where dQ"^ is the line element of the 2-sphere and all the four functions A^, ip, a, f3 depend on 
t and r. We shall work in isotropic coordinates [jl^ (a^ = ip'^). Accordingly we shall have 
only one conjugate variable to '0 by a diagonal gauge in which the shift vector [3 is set to 
zero: in this way the extrinsic curvature will be determined only by its trace ^7^. Therefore 
the metric we shall use is the following 

ds^ = -N\t, r)de + ^^(t, r) \dr^ + r^dfi^j _ (3) 



The equation of motion for a scalar field with a generic potential V{(f)) 



dV 

- V''d,<l), + TT- = (4) 



are written in a Hamiltonian form for the metric (^ 



. TT,: 



^# (^) 



.dV 

r 



1 {rVN^,^ - W^ (6) 



where 0, ttj are the pair of conjugate variables for the scalar field, a dot denotes the derivative 
with respect to t whereas a prime denotes the derivative with respect to r. 
The Hamiltonian equations for the metric are the following: 

i---^K (7) 

where K is the trace of the extrinsic curvature and E is the total Hamiltonian for the scalar 
fields 

The equation for the lapse function A^ does not involve time derivatives (since A^ is not a 
dynamical degree of freedom) : 

The energy and momentum constraints are 



i'"«--^vV 
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(11) 



lK' = y^K^^ (12) 

3 Y ^' 

where ^^^R is the curvature of spatial sections. 

III. QUASILINEAR REGIME 

To describe analytically the phase that follows the linear stage, we employ now a quasi- 
linear approximation, which generalizes the method developed in 0,0 for rigid spacetimes. 
In this approximation we consider terms that are linear in the fluctuations of the inflaton 
and of the metric, but keep quadratic terms in x fluctuations. In the absence of an initial 



homogeneous value of X; these quadratic terms represent the leading effects of x o^^ the 
infiaton and the metric functions. 

The metric at the first order has the same general form as in the linear theory (reviewed 
in Ref. [Q)- Including only scalar perturbations but keeping (for now) three-dimensional 
coordinates, we write it in the longitudinal gauge as 

ds2 = a^(t)|-[l-2$(t,x)]rft2 + [l + 2^(t,x)]rfx2} ; (13) 

t now is the conformal time. An important difference with the linear theory, however, is that 
even for minimally coupled scalars, there is an anisotropic stress, induced by a nonlinear 
contribution of x- As a result, $ 7^ \E'. 

For definiteness, from now on we consider a model with two fields: the infiaton (p and 
the second field x whose fluctuations will drive fluctuations of and the metric. Both fields 
are minimally coupled to gravity, and their potential is 

n</',x) = ^A0^ + ^/0V + ^AV. (14) 

Our method applies equally well to single field models, in which case the role of x will be 
played by short-scale components of the infiaton itself. 

To the linear order in metric perturbations, equations for the fields and x can be 
obtained without restriction to spherical symmetry. On the other hand, equations for the 
metric functions we present below are valid only under that restriction. 

Expanding the equation of motion for to the first order in $, \1', x^, and the spatial 
derivatives of we obtain 

+ 27^0 + V^0 + AaV^ + ^^aVx^ = -<P{^ + 3^) + 2Xa'^(f)^^ , (15) 

where TC = d/a; dots denote derivatives with respect to t. We can separate the field into 
its volume average, 0o, and a fluctuation: 

0(t,x) = 0o(t) + 50(t,x). (16) 

The volume averages of $ and "^ are made zero through definitions of the scale factor a and 
the conformal time. We also assume that the volume average of the field x is initially zero. 
Substituting (|T6|) into (|1^), expanding to the first order, and averaging over the volume 
gives 

00 + 27^00 + Xa^^l + g^a^^o¥ = ; (17) 



the overline denotes a volume average. Subtracting {\n\) from (|T5|) , we obtain an equation 
for the fluctuation: 

50 + 2n6^ + [-V^ + 3Xa^(f)l\ 6(f) + g^a^Mx^ - F) = -M^ + 3^) + 2Xa^(f)l^ . (18) 

Similarly, for x, to the first order in $, \E', x^, and 50, we obtain 

X + 2nx+ \-V' + g\^<j)l\ X = -2g^a^(t>oXH - x{^ + 3^) . (19) 



In this paper, we use quasilinear theory only for the earhest times, when the right-hand 
of ([T9|) can be neglected. In this case, we obtain an equation for x which is unaffected by 
metric perturbations: 

X + mx+ [- V^ + g^a^(t>'S\ X = . (20) 

In the full 3d case, the nonlinearities of x would act as a source not only for scalar metric 
perturbations, but also for vector and tensor ones. Equations given below are restricted to 
spherical symmetry, in which case only scalar perturbations are present. 

In the quasilinear approximation, the energy constraint (^) reads 

- v^m + ?,n{^ + n^) = — {6E4, + 6E^) (21) 

where the fluctuations of energy densities of the fields are defined by 5Ei = Ei — Ef, for the 
model at hand, to the required order, 

6E^ = <Po6<P + 0^$ + a^X(t)l64> , (22) 

i^x = Y + ^ + J-'<Plx' . (23) 



The momentum constraint (|T2D reads 

(vi/ + 7^$)' = -^(0o50' + xx')- (24) 

We do not present a quasilinear form for Eq. (P). The quasilinear form for Eq. (|10]) is: 

.((i^)'..^,^ (25) 

The right-hand side is the anisotropic stress. 

In general, an anisotropic stress is generated by a non-vanishing traceless part of the 
spatial energy-momentum tensor. In our case, it is due to the traceless second-order contri- 
bution of the field x'- 

n} = d'xdjX - l^pxdiX ■ (26) 

For small scalar perturbations, the corresponding Einstein equation takes 

{D'Dj - -6)0^) [$ - ^] = ^2 [ml , (27) 
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where Di is the covariant derivative with respect to the spatial metric jij, related to gij via 
dij = ip^lij, D^ = D^Di, and the subscript S denotes the scalar part of 11*. The scalar part 
is automatically selected when we make restriction to spherical symmetry, and the result 
coincides with Eq. (|25|) . 

We see that nonlinearities in x can potentially cause two types of effects. First, the 
rescattering term — the last term on the right-hand side of ([T8| ) — will lead to a growth of 50; 
this is the effect represented by the diagram of Fig. 1. Second, the right-hand side of Eq. 
(^) directly produces $ — \l/ 7^ 0. Either way, the resulting perturbations are non-Gaussian. 
In the remaining sections, we will study these two effects in turn. Note that neither of these 
effects was taken into account in Ref. |jl9|] where only the direct second-order contribution 
of X to curvature perturbation was considered (in a different model) . 
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IV. INFLATON FLUCTUATIONS 



The two effects described at the end of the previous section are more easily disentangled 
if we use, instead of 6(j), the gauge-invariant variable 



-<j)'^/n 



(28) 



This variable was introduced by Mukhanov to characterize perturbations produced during 
inflation [20|; it was applied to studies of reheating at the linear level in Refs. mT|l2-14|.P] 
From now on we omit the zero subscript on the homogeneous component of 0, hoping that 
this will not cause any confusion. 

Using Eq. (|TBp and the quasilinear equations for the metric, we obtain the following 
equation for Q^: 



Q^ + 2nQ^ + [-V' + 3AaV' + ^grav] Q'P + ^'a'0^x' = F{x) , 
where 6x^ = X^ ~ X^-, the gravitational contribution djQj^Av ^^ ^^^ frequency is 



(29) 



"^GRAV 



.«^t. 



K 



27e 



(30) 



and V is the potential of 0; in our case, V^^ = X(f)^. This potential does not include a Hartree 
correction, cf. Eq. (0). Indeed, including such a correction here would take us outside of 
the quasilinear approximation. Thus, (|30D is the standard expression of the linear theory 
8|JT^. The force term F{x) depends only on x'- 



Fix) = f^' 




(?V0%^ + — V'($-vl/) 



3k^ 



K 



n^ 



z. 



X ' 



(31) 



where Z^ is the contribution of x to the curvature perturbation, see eqs. (|§), ( ^Tj ) be- 
low. The term proportional to V^($ — \1') depends only on x by virtue of the quasilinear 
approximation, see Eq. (p5|). 

Because Eq. (|29D, unlike Eq. (p!8D, depends on the gravitational potentials only through 
V^($ — ^), the most long-range (secular) contributions to $ and \1' do not enter, and all 
interactions generated by gravity in Eq. (^) are suppressed by at least one power of Mpi. 
(This suppression is in parallel with the result of Ref. ||13| for the linear case.) Neglecting 
all such terms and going over to the variable X = aQ^, we obtain 



X + 



-V^ + 3\a^d)^)X 



g'aMx' - X') 



(32) 



^This variable looks as a combination of metric and field fluctuations in the longitudinal gauge, 
while in the uniform curvature gauge its interpretation is more transparent: it is simply the field 
fluctuation, which indeed is the true dynamical variable of the problem (see for example p2|). 



where we have neglected the term a/a in the frequency because a radiation dominated 
universe (a ~ 0) is a good approximation during reheating in this model. Eq. (|32D shows 
that, in the quasilinear theory, the growth of X is driven primarily by rescattering due to 
the local coupling and is essentially unaffected by gravity. 

By using the spherically symmetric equations for gravity in deriving Eq. (p2D, we have 
neglected vector and tensor perturbations. It seems plausible, however, that contributions 
of these types of perturbations to the evolution of Q^ are suppressed in the same way as 
the contribution of scalar perturbations. We therefore expect that Eq. (^) applies, to the 
leading order in Mpj^, in the full three dimensions. Accordingly, we write it in the Fourier 
space as 

Xk + (A;2 + ^\a'(f')X^ = -g^a^ J d'pxpit')xu+pit') . (33) 

The solution to Eq. ( ^31) can be written in terms of the Green function: 

Xk(t) = Xk,iNFL(t) - g' J^ dt'Gkit, t')a\t')4>{t') J d^px;{t')xi.+p{t') , (34) 

where Xk,iNFL(^) is the solution to the linear problem, corresponding to perturbations pro- 
duced during inflation and the crossover to reheating. The second, forced, term in Eq. (|3^ 
starts from zero at t = but then grows, initially as exp(2/iinax^) where //max is the maximal 
Floquet exponent for x- There is no additional decay in time, since the factor a^ cancels 
the damping envelope ~ 1/a of the fields 0, x- 
The Green function can be written as 

(^k[t,t) = — , (35j 

Wk 

where XI-I2 are two independent solutions to the homogeneous equation 

Xi%it) + (P + 3Xa'^')xi% = , (36) 

and Wk is their (time independent) Wronskian. The homogeneous solutions are known in 
imphcit form |^ (see also p4|): 



Xn,M^^m)\e.v(±'^lj^). (37) 



where Mj^{t) is defined in Eq. (37) of Ref. [Q, and the Wronskian W^ is 



iy, = ^fc^(x4-^)(3-x4) (38) 

with X = k/{V\(l){0)). 

The Wronskian Wk is real in the single resonance band 3/2 < x^ < ^/S and imagi- 
nary everywhere else. When Wk is imaginary, Wk = i\Wk\, the solutions ( pTf ) are complex 
conjugate, and the Green function becomes 



Gkit.t') = J \Mk(t)Mk(t')\ sin ^—^ / . (39) 

At small fc, \Wk\ ~ ^; while Mi^[t) is finite and nonzero, so G^ is regular there. Thus, if 



the field x has not been amplified during inflation, the integral over p in Eq. (34) remains 
regular at A; = as long as the quasilinear approximation applies. 

Using Eq. (|3^) , we can estimate the power spectrum of the forced (non-Gaussian) part 
of Q^ as 




Pq,ng(A;) ~ f7 / ^>^xb)^x(IP + k|) , (40) 



where P^ is the power spectrum of x- We should now distinguish two cases. 

(i) If the field x was light during inflation, g^ < A, the spectrum of x produced at that 
time is infrared singular, which by Eq. ( |40| ) gives rise to an infrared-singular spectrum of Q,^. 
This will result in non-Gaussian adiabatic perturbations of the metric at scales of relevance 
to structure formation, in addition to the previously considered [^ perturbations due to a 
direct quadratic contribution of x to the energy density. 

For 1 < g^ jX < 3, when the infrared part of the spectrum is amplified by parametric 
resonance, this effect can be quite large. It is the largest for g'^/X = 2, when the resonance 
peak is centered at A; = 0. In this case, assuming for simplicity the exact scale-invariant 
spectra for both and x at the end of inflation, we can estimate P^ as Px{k) ~ (x^)fuii/^^) 
where (x^)fuii denotes the full 3d variance. The resulting Pq,ng should be compared to 
the inflation-produced Gaussian spectrum Pqjnfl ~ -^infl/^^; where -^infl is the Hubble 
parameter during inflation. Using the maximal value of the variance (x^)fuii ~ Xcj)^ / g"^ [0)11; 
we obtain 



2 - „2 -lOV (41) 



-Pq,INFL \M ) -f^INFL -f^INFL 

(A similar, but Gaussian, effect is obtained in linear theory in certain models with a nonva 



nishing homogeneous component of x [13,1 



(ii) If the field x '^as heavy during inflation, it arrives at reheating with a constant 
(white-noise) spectrum at small fc, and by Eq. ( ^OD the spectrum of Q^ is also constant 
("non-Gaussian white noise"). Indeed, assuming that the resonance peak is located at 
k = /c*, we estimate that P^ ~ {x'^){uii/kl, where (x^)fuii denotes the full 3d variance. Then, 

Using again the maximal value (x^)fuii ~ Xcp'^/g^, we obtain 

Pq,ng ~ ^Vkl . (43) 

This estimate can be converted into an estimate for curvature perturbations, as measured 
by the parameter (, 



where 



[l + w)C='^(4//n + '^)+{l + w)<i' , (44) 



l + -=3^. (46) 



The term 1 + w in front of ( is needed in order to avoid singularity at </> = |[r2| , p^ . Since 
1 + w just oscillates between finite values, a growth of the left-hand side of (0) can only 
come from a growth of (. 

The momentum constraint (p^, can be generalized to the full three dimensions; the 

is 



resulting scalar part, separated along the lines of Ref. p6 

^k + H$k = -y (c/'^^k + ^ / dM^ ■ p)Xp+kX;) = - Y#</>k - ^xk • (46) 

Using (|46|), we can express C as 

(l + z/;)Ck = -^0g^k- J^^xk- (47) 



At small fc, we estimate the integral in (46), in case (ii), as 



|d='p(k-p)xp+kX;~C)(fc'), (48) 

so that the direct contribution of x i^ this case is independent of k in the small- A; limit. 
The same conclusion was reached previously in Ref. [|19|, albeit for a different model and 
by a different method. However, comparing (|^) and (|40|), we see that at 5^^ :^ A the 
contribution from Q^, which was not considered in Ref. |jl9|, is parametrically larger than 
the direct contribution of x- 

Even that larger contribution is still quite small. Using (^31) , we obtain 

SO that the perturbation of ( on scales larger than the Hubble scale is 

where we have used 0/Mpi ~ 3 x 10"'^ and Ti/k^, ~ 6 x 10"'^. We doubt that such a small 
perturbation can cause a gravitational instability, such as one observed numerically in Ref. 
||16|| . Indeed, in the following sections we will identify a much stronger infrared effect, due 
to anisotropic stress. 

The form (^) of the Green function leaves open the possibility of 6(j) cascading towards 
the infrared. Indeed, if the integral in the argument of sine increases with time, the region of 
k where the regular small-fc limit applies will be pushed further and further in the infrared. 
However, the corresponding growth of the spatial scale of the perturbations can only occur 
at the relevant excitation speed, such as the speed of sound, and so cannot overtake the 
cosmological horizon (which recedes at the speed of light). 
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V. INITIALIZATION OF A RADIAL FIELD 



It is clear from Eq. ( pSj ) that to understand the effect of the anisotropic stress, we need 
to understand the spatial profile of the field x- If is useful to look first at the initial state, 
because the profile of x remains qualitatively the same throughout the quasilinear stage of 
the evolution. 

A spherically-symmetric (also called radial) field x is initialized according to the formula[| 

rx{r)= AY^ SkkP^^'^ik) sin kr , (51) 

k 

where k are the lattice momenta, quantized in units of tt/R, Sk are Gaussian random numbers 
with zero average and unit dispersion, and P^. is the power spectrum. The constant coefficient 
A is to be chosen upon deciding what part of the full three-dimensional (3d) variance the 
radial field will account for. The latter question arises because in the full 3d case there surely 
will be other, nonradial, modes, which will carry a part (in fact, most) of the variance. 

Consider the average of the field defined in (^) over the distribution of the random 
numbers s^: 

{X^)(r) = — - / dk^P^ik) sin^ kr , (52) 

where we have replaced the sum with an integral, and fcmin and fcmax are the minimal and 
maximal lattice momenta. For a field that was heavy during inflation, we use Px{k) corre- 
sponding to the adiabatic vacuum, with an additional ultraviolet cutoff: 

P^ik) = {2n)-''{2uJk)-' exp{-e/kl) . (53) 

where Uk is the instantaneous frequency at t = 0. For this power spectrum, the average ( |5^ 
has the asymptotic form 

{x'){r)^^^JPAk)k'dk (54) 

at large r. On the other hand, the full 3d variance is independent of the coordinates and 
equals 

{X^)iuii = ^7r J Px{k)edk. (55) 

The field x is smooth, almost constant, over a patch of radius r ^ vt/Zcq, where fco is 
the ultraviolet cutoff from (^5]). A smooth field is predominantly an s-wave, so for such r 
we expect the average (|52D to approximately equal the full 3d variance (|55|). In a sense, the 
scale 



^ In numerical work, after the field x has been initialized according to (O) we additionally subtract 
away its homogeneous component. This does not matter for the present argument. 



11 



4 = vr/fco (56) 



can be viewed as the initial correlation length of the field. The asymptotic expression (54) 
applies at r ^ Eq, but as r approaches io, we expect {x'^){r) to approach (x^)fuii- Accordingly, 
we adopt the following estimate for the magnitude of (x^)(^) at r ^ Eq: 

ix'Kr) = {X%unil/r' . (57) 

This corresponds to 

A^ = 8itX/R ■ (58) 

A number of order one could be inserted here, but it has practically no effect on the final 
results. The important point is that, at a fixed r, Eq. (|^) gives the size of fiuctuations that 
is independent of the spatial volume. Equivalently, the volume average 

yj\'dr{x'){r) = ^-^{x'Uu (59) 

decreases with R as l/i?^. 

The 1/r^ dependence of the variance of a radial field can be simply interpreted as follows. 
For an observer at the origin, a smooth patch of size io at distance r ^ Iq has angular size 
io/r. Thus, expansion of a 3d field with correlation length io will require of order {r /CqY 
different spherical harmonics, with amplitudes that are roughly of the same order. The 
variance in each of these partial waves, specifically the s-wave, will then account for {(-o/rY 
of the full 3d variance. Similar considerations apply to perturbations of the metric. 

In Ref. |T^, in preparation of the initial state, we in effect assumed that the volume 
averages of (x^) and of the variances of the metric functions, e.g. the left-hand side of (^9|), 
equaled the full 3d variances. As we now see, that was an overestimate of the size of initial 
fiuctuations. However, as typical in problems where amplification of fiuctuations is cut off 
by nonlinear effects, the initial size is relatively unimportant: repeating the run of Ref. W^ 



with initial data of the correct size, we find that the strong gravitational effects, which we 
interpreted as formation of a black hole, are merely postponed from t ~ 74 to t ~ 105 (in 



the units of Ref. \lt 



VI. EFFECT OF THE ANISOTROPIC STRESS 

We are now in position to estimate metric perturbations caused by the anisotropic stress. 
Consider Eq. (p5D of the quasilinear approximation. As we have seen in the previous section, 
the initial profile of x is a superposition of spherical waves and has a 1/r envelope. This 
form is preserved during the quasilinear stage for as long as Eq. (pOD applies. However, the 
correlation length of x is now determined not by the ultraviolet cutoff, but rather by the 
typical resonant momentum, k,,: 

I = n/K ■ (60) 

Accordingly, instead of (^) we now have, at r ^ £, 
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{x'){r,t) = {x'hMty'{t)y. (61) 

Averaging Eq. ( pSf ) over the distribution of the random initial data and using Eq. (|6l|) , we 
obtain 



m-{^)y 



c^'^^ , (62) 



where c > is a weakly time-dependent constant of order one. Two powers of i in Eq. (|6TD 
eat away two spatial derivatives in Eq. (^). 

Eq. ( |62D is easily integrated. Using the boundary conditions ($ — '^){R) = 0, ($ — 
^)'(i?) = 0, we obtain^ 

($) - (^) = Ic^'ix'Ui {in J + l"-^} ■ (63) 

Thus, the average gravitational potentials are more long-ranged than the average of x^ itself, 
Eq. (^. 

Combining Eq. (|63D with the momentum constraint (|5), we can obtain estimates for 
($) and (\&) separately, at least for r > i. We will need to integrate Eq. (^) so as to 
remove the spatial derivative. The second term on the right hand side can only lead to a 
short-range, i.e. 1/r^, contribution, while the first term leads to a contribution that is either 
short-range or suppressed by several powers of Mpi. So, for our present purposes we can 
take 

(^) + H($) ^ . (64) 

If, as a result of parametric resonance, the field x grows with a characteristic exponent /x, 
the potentials $ and \E', as (^) and (|6^) show, grow with characteristic exponent 2/i. Then, 
from (|6^) we see that 

(^)^^($)<($), (65) 

/imo 

where mo ~ vA</)(0) is the frequency (in conformal time) of the infiaton's oscillations. 
Therefore, 

W^i^lnf, (66, 

and then (|65|) is an estimate for (\E'). We see that, up to a weak logarithmic factor, the 
average potentials are independent of R. We stress, however, that these estimates apply 
only in quasilinear theory. 



Using, instead of ($ — ^){R) = 0, the condition of vanishing volume average of <I> — ^, as in 



Sect. [D, will only change ( p3D by an unimportant constant. 
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The approximate volume independence of ($) and (\l/) during the resonant growth of x is 
confirmed by numerical integrations. Integrations were done using the algorithm described 
in Ref. |T^ and are fully nonlinear. For example, Fig. |^ shows Sn'ip/'ipiR) for g'^/X = 110 



and A' = 0. Here ijj{R) is the conformal factor at the outer end of the grid, and Sn'ip is 
defined by 



rn 



[Snipf = 47r /" dkk^P^ik) , (67) 

where P^{k) is the power spectrum oi ip — ip{R), and k^[^ is the minimal lattice momentum. 
Fourier components oi ip — ip{R), required for calculation of the power spectrum, are deter- 
mined via a sine transform. During the resonant growth of x? the metric functions $ and 
\1' rapidly become smooth and can therefore be regarded as self-averaging. For this reason, 
the quantity plotted in Fig. |^ should give a good idea of the typical values of (\1/). 

Numerical values of ($) and (\1/) predicted by the estimates ( |66|) and (|^) are quite small. 
The usual estimate [0,|[ for the maximum value of (x^)fuii is 



A0' A lO-'M] 



2/1 #2 
PI 



max(x')fuii ^ ^ X , /' , (68) 

g^ g^ In A 

so for A = 10^^'^ and g^/X ~ 100 we obtain 

(*)<($) ~ 10-Mn(i?/£) . (69) 



From Fig. g, we see that the lattice result by far exceeds the estimate (|69|), even during the 
stage of exponential growth, which one might want to identify with the quasilinear stage. 
The reason is that, because of the 1/r^ form of (x^), the quasilinear approximation, strictly 
speaking, breaks down quite early on — when (x^) becomes of order ( |68|) at r ~ £. The 
growth of various quantities can still continue, however, because in most of the space (x^) 
is still much smaller than that. (We discuss this point further in the next section.) 

During the entire exponential growth, the two curves corresponding to different box sizes 
coincide with excellent precision, in agreement with the approximate volume independence 
of (\1/) derived analytically. The time when the curves part corresponds to breakdown of the 
quasilinear theory in most of the space. At that time, volume dependence of the results is 
not surprising: nonlinear gravitational phenomena, such as gravitational collapse, depend 
on the total amount of matter available and will be more pronounced in a larger volume. 



VII. FURTHER NONLINEAR EFFECTS 

The driving force due to local scattering in Eq. (|18]) can be estimated as 

g'<P{x')-g'<l>f{x'huiiy. (70) 

Fluctuation of (p caused by this driving force at r ~ ^ is 



^ - p^' 
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This becomes of order when (x^) reaches the value (0). After that (x^) does not grow 
any further at r ~ £ but continues to grow at r 3> ^. 

From our lattice data, it appears that this later growth is a combination of the usual res- 
onant growth and rescattering of fluctuations from smaller r, where they are already large.. 
As a result, the region where (x^) is large propagates to larger r, so that the quasilinear 
form (|6TD no longer applies. From Eq. (^51) , we see that the growth of this region leads 
to a further growth of metric perturbations. If the total volume is large enough, these can 
eventually become sufficiently large to cause a gravitational instability. 

The time when the resonant growth of x has ended in most of the space corresponds to a 
short plateau around t = 250 in Fig. |^. The subsequent growth of Sn'ip has a different origin: 
we interpret it is a result of an additional production of 6(f) by another, weaker resonance, 
due to the A0^/4 interaction. During this further growth, we observe formation of a black 
hole — the lapse A^, as a function of r, crosses zero.^ 

The presence of the second stage of growth in Fig. ^ leads to an important conclusion. If 
some decay channel becomes blocked while the infiaton's homogeneous component still has 
not decayed completely, the infiaton will "wait" until a new channel opens up. At least in 
our spherically-symmetric geometry, and in the model considered here, its decay via the new 
channel still causes significant metric perturbations. This conclusion holds also in the case 
when the decay of the infiaton into x is blocked by a A'x^/4 self-interaction (cf. Ref. 0): we 
show the result for g'^/X = 8 (a stronger resonance) and A'/A = 1000 by short dashes in Fig. 
^. In this case, the decay into x is cut off at much smaller values of metric perturbations, 
and the interpretation of the second stage of growth as the A0^ resonance is unambiguous. 

VIII. CONCLUSION 

In this paper we have studied the influence of nonlinear effects, arising during decay of 
the infiaton, on metric perturbations. Although we restricted our attention to a quartic 
potential for the infiaton, results for other models of chaotic infiation can be obtained by 
similar methods, and we believe they will be similar in nature. 

Early stages of the infiaton decay can be studied analytically, using the quasilinear 
approximation, at least under the restriction to spherical symmetry (a similar approximation 
may exist also in the full 3d case). The quasilinear approximation retains quadratic terms in 
the decay products; these are the leading terms affecting metric perturbations after inflation 
in the case when flelds interacting with the inflaton do not have sizeable homogeneous 
components. 



^ One may, and some physicists did, debate whether this is an adequate definition of a black 
hole. Our main argument remains the same as in Ref. [16|: for the metric (^, a sphere on which 



A^ = cannot be crossed by a classical particle, and is therefore a horizon. On the other hand, 
quantum particles, or classical waves, can transport energy across the horizon (similarly to Hawking 
radiation); we suppose this transport is the cause of growth or evaporation of black holes in our 
simulations. 
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Using the quasilinear approximation, we have identified two effects that could potentially 
produce large-scale metric perturbations after the end of infiation. We have found that 
the first of these effects — scattering of the second field x off the infiaton's homogeneous 
component — does not cause significant super-Hubble perturbations in the quasilinear regime. 
However, the second type of effect — metric perturbations caused by the anisotropic stress — 
has been found to be important. 

It is interesting to contemplate how this latter effect is able to avoid the restriction on 
super-Hubble physics imposed by causality. To this end, note that a spherical wave has a 
large degree of angular coherence: even for a wave with momentum much larger than the 
Hubble parameter, angular coherence extends over scales that are not in causal contact. 
Parametric amplification of such fiuctuations is not restricted by causality (just as it is not 
restricted during infiation). 

Persistence of angular coherence of fiuctuations in the full three dimensions corresponds 
to a large mean- free path. As we have discussed in Sect. |I|, this is indeed likely to be the case 
at the early stages of reheating in weakly-coupled models. However, because restriction to 
spherical symmetry explicitly breaks translational invariance, it is not immediately obvious 
how our results can be transplanted to the full 3d case. Indeed, the r dependence of the 
ensemble averages ($) and (\1') (cf. (|53D) is a direct consequence of the translational non- 
invariance and is not expected to persist in three dimensions. Instead, the 3d picture that 
we have in mind is a sort of "metric turbulence" , in which a (nearly) spherically-symmetric 
fiuctuation in a sphere of radius -Ri sits on top of a fiuctuation in a sphere of a larger radius 
i?2- Then, instead of ($) and (\E') we should consider the total perturbation from all scales 
shorter than some i?, e.g. 

5fi$ = ([$(x)-$(x + R)]2)V2^ (72) 

so for example the quasilinear result (p^) will correspond to 



(with an approximately scale-invariant 3d spectrum Pq,{k) oc A;^'^lnA;). 

Perturbation (|73D is by itself quite small. We have seen, however, that in spherical 
symmetry a much larger perturbation is possible, due to spreading out of the region where 
fiuctuations of x ^"^^ large. In Sect. |VII| , we have argued that this is indeed the mechanism 
of gravitational instability and formation of black holes observed in our lattice studies. 
Because this mechanism involves already large field fiuctuations, we cannot exclude the 
possibility that in the full three dimensions scattering of these large fiuctuations destroys 
angular coherence, by breaking the s-wave into higher partial waves, and the formation of 
black holes is prevented. However, if black holes do form in the full 3d by the mechanism 
identified here, the question remains what is the spectrum of metric perturbations after 
these black holes evaporate. 

This work was supported in part by the U.S. Department of Energy through Grant 
DE-FG02-91ER40681 (Task B). 
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FIGURES 





(a) 



(b) 



FIG. 1. Diagrams representing effects of the products x of the inflaton decay on metric pertur- 
bations: (a) through inflaton fluctuations and (b) directly. 
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FIG. 2. Relative variation of the conformal factor over super-Hubble scales, as a function of 
conformal time (which is in units of l/vA0(O)) for g"^ /X = 110, X' = 0, R = SOOtt, number of 
grid points n = 8192 (solid line); g'^/X = 110, X' = 0, R = IOOOtt, n = 16384 (long dashes); and 
c,2/A = 8, X'/X = 1000, R = SOOtt, n = 8192 (short dashes). 



